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(54) Frequency shift examination in a CDMA receiver 


(57) The invention relates to a method for acquiring 
a receiver (1) into a code modulated spread spectrum 
signal received by the receiver (1 ). In the method at least 
one replica code (r) is used which corresponds to a code 
used v/ith the modulation having a pre-determined 
number of chips, and an examination phase is per- 
formed, in which a frequency shift of the signal within a 
selected frequency area is examined, and a code phase 
of the code used with the modulation is examined. In the 
method the frequency shift examination is divided into 


a first estimation phase and a second estimation phase, 
wherein in the first estimation phase the selected fre- 
quency area is divided into a first set of frequencies, and 
in the second estimation phase a second set of frequen- 
cies is examined nearby each frequency of the first set 
of frequencies. Further, a comparison for frequencies of 
the second set of frequencies is performed using the re- 
ceived signal and the replica code (r). and the results of 
the comparison are used to estimate the correct fre- 
quency shift. 
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Description 

[0001] The present invention relates to a method for performing acquisition in a receiver as set forth in the preamble 
of the appended claim. 1. The present invention relates further to a location system as set forth in the preamble of the 

5 appended claim 8. The present invention relates further to a receiver as set forth in the preamble of the appended 
claim 15. The present invention relates also to an electronic device as set forth In the preamble of the appended claim 22. 
[0002] One known positioning system is the GPS system (Global Positioning System) which comprises more than 
30 satellites, of which at least 4, in sorrie situations even 12, are simultaneously within the sight of a receiver. These 
salelliles transmit information about satellites, e.g. Ephemeris data of the satellite as well as data on the time of the 

10 satellite. The receiver used for positioning infers its position normally in such a way that the receiver calculates the 
time of propagation of a signal transmitted simultaneously from several satellites in the positioning system to the re- 
ceiver For the positioning, the receiver must typically receive the signal of at least four satellites within its sight, in 
order to calculate the position. 

[0003] Each satellite of the GPS system transmits a so-called LI signal at a carrier frequency of 1575.42 MHz. This 

'5 frequency is also indicated with 154to. where f^ = 10.23 MHz. Furthermore, the satellites transmit an L2 signal at a 
earner frequency of 1227.6 MHz, i.e. 120fo. In the satellite, the modulation of these signals is performed with at least 
one pseudo sequence. This pseudo sequence is different for each satellite. As a result of the modulation, a code- 
modulated wideband signal is generated. The modulation technique used makes it possible in the receiver to separate 
the signals transmitted from different satellites, although the carrier frequencies used in the transmission are substan- 

20 tially the same. This modulation technique is called code division multiple access (COMA). In each satellite, for mod- 
ulating the LI signal, the pseudo sequence used is e.g. a so-called C/A code (Coarse/Acquisition code), which is the 
Gold code. Each GPS satellite transmits a signal by using an individual C/A code. The codes are formed as a modulo- 
2 sum of two 1023-bit binary sequences. The first binary sequence G1 is formed with a polynome X^o + + 1, and 
the second binary sequence G2 is fomned by delaying the polynome X^^ + X® -t- X^ + + + X^ + 1 in such a way 

?5 that the delay is different for each satellite. This arrangement makes it possible that different C/A codes can be produced 
with an identical code generator. The C/A codes are thus binary codes whose chipping rate in the GPS system is 1 .023 
MHz. The C/A code comprises 1023 chips, wherein the iteration time of the code is 1 ms. This repeating sequence 
can also be called as an epoch. The carrier of the LI signal is further modulated with navigation information at a bit 
rate ofSObit/s. The navigation information comprises information about the health of the satellite, its orbit, time data, etc. 

30 [0004] During their operation, the satellites monitor the condition of their equipment. The satellites may use for ex- 
ample so-called watch-dog operations to detect and report possible faults in the equipment The errors and malfunctions 
can be instantaneous or longer lasting. On the basis of the health data, some of the faults can possibly be compensated 
for, or the information transmitted by the malfunctioning satellite can be totally disregarded. Furthermore, in a situation 
in which the signal of more than four satellites can be received, different satellites can be weighted differently on the 

35 basis of the health data. Thus, it is possible to minimize the effect of errors possibly caused by satellites which seem 
unreliable in the measurements. 

[0005] To detect the signals of the satellites and to identify the satellites, the receiver must perform synchronization, 
whereby the receiver searches for the signal of each satellite at the time (acquisition) and attempts to be synchronized 
and locked to this signal (tracking) so that the data transmitted with the signal can be received and demodulated. 

40 [0006] The positioning receiver must perform the synchronization e.g. when the receiver is turned on and also in a 
situation in which the receiver has not been capable of receiving the signal of any satellite for a long time. Such a 
situation can easily occur e.g. in portable devices, because the device is moving and the antenna of the device Is not 
always in an optimal position in relation to the satellites, which impairs the strength of the signal coming to the receiver 
AlsO; in urban areas, buildings affect the signal to be received, and furthermore, so-called multipath propagation can 

45 occur, wherein the transmitted signal comes to the receiver along different paths, e.g. directly from the satellite (line- 
of-sight) and also reflected from buildings. This multipath propagation causes that the same signal is received as 
several signals with different phases. 

[0007] Both the movement of the receiver and the movement of the satellites in general will continuously change the 
strength, phase, number, etc. of the signals entering the receiver The movement of the satellites and Inaccurate local 

50 time reference of the receiver also cause a shift of the carrier frequency, also called as a Doppler shift. Thus, changes 
are also caused in the signal received by the receiver. For this reason, the receiver must try. after the synchronization, 
to continuously maintain locking to the signal of each satellite from which information is received e.g. for positioning. 
In the receiver, the code phase ts calculated very frequently, and the oscillator is adjusted, if necessary, in such a way 
that the receiver remains synchronized. 

55 [0008] The positioning arrangement has two primary functions: 

1 . to calculate the pseudo distance between the receiver and the different GPS satellites, and 

2. to determine the position of the receiver by utilizing the cateulated pseudo distances and the position data of the 
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satellites. The position data of the satellites at each time can be calculated on the basis of the Ephemeiis and 
time correction data received from the satellites. 

[0009] The distances to the satellites are called pseudo distances, because the time Is not accurately known in the 
5 receiver. Thus, determinations of the position and the time are repeated, until a sufficient accuracy has been reached 
with respect to the time and the position. Because time is not known with absolute precision, the position and the time 
must be found out e.g. by linearizing the set of equations for each new iteration. The pseudo distance can be calculated 
by measuring the mutual, virtual propagation delays of the signals of the different satellites. 

[0010] Almost all known GPS receivers utilize correlation methods for calculating the distances. In a positioning 
10 receiver, pseudo random sequences of different satellites are stored or generated locally, A received signal is subjected 
to conversion to an intermediate frequency (down conversion) and then to the baseband, after which the receiver 
multiplies the received signal with the stored pseudo random sequence. The signal obtained as a result of the multi- 
plication is integrated or lowpass filtered, wherein the result is data about whether the received signal contained a 
signal transmitted by a satellite. The multiplication is iterated in the receiver so that each time, the phase of the pseudo 
'5 random sequence stored in the receiver is shifted. Consequently, this means cross -correlation between the received 
signal and the pseudo random sequence generated/stored in the receiver. The correct code phase is inferred from the 
correlation result preferably so that when the correlation result Is the greatest, the correct code phase has been found. 
Thus, the receiver is correctly synchronized with the received signal. 

[001 1] This code synchronization process should be performed at different possible Doppler frequencies, otherwise 
20 even for correct code phase the correlation could be small if the received signal is modulated due to Doppler shift In 
frequency. 

[0012] After the synchronization with the code, the next step is the phase locking. The congelation result also indicates 
the information transmitted in the GPS signal; that is, it is a demodulated signal. 

[0013] The above-mentioned synchronization and frequency tuning process must be iterated for each signal of a 
25 satellite which is received In the receiver Consequently, this process takes a lot of time particularly in a situation in 
which the signals to be received are weak. To speed up this process, some receivers of prior art use several correlators, 
wherein it is possible to search for correlation peaks simultaneously. However, in practical applications it is not feasible 
to speed up the acquisition process by increasing the number of correlators because the number of correlators can 
not be increased unlimitedly due to e.g, the power consumption, the size and the cost of the receiver. 
30 [0014] In some prior art receivers FFT technique have been used in connection with the correlators to define the 
Doppler shift of the received GPS-signal. In these receivers a con-elation is used to decrease the bandwidth of the 
received signal. Then this narrow-band signal is analyzed with FFT algorithms to estimate the Doppler shift in carrier 
frequency. 

[0015] It is an aim of the present invention to provide a receiver in which an attempt is made to estimate the Doppler 

35 frequency and code phase shifts for acquisition to the received signal. The invention is suitable for use particularly in 
positioning receivers, but also in other receivers, preferably CDMA receivers, where the receiver must be synchronized 
and locked with a spread spectrum signal. The invention is based on the idea that the search in Doppler frequency 
dimension is performed for a group of nearby frequencies at once along with the search in code phase dimension. 
Thus the search in Doppler frequency is performed essentially in two phases, wherein in one phase a coarse Doppler 

'*o frequency is assumed and then a fine grid of Doppler frequencies is analyzed around that coarse frequency. At the 
same time the search in code phase direction could be performed by many known techniques implementing circular 
correlation. The preferred approach described In this invention assumes the correlation is performed in frequency 
domain using D FT (FFT). In the present invention the Doppler frequency search is actually perfor med by compensating 
t he^os^ibjej reguenc^ s hl ft, then correlating with locafireplicac ode for different codeptras^sand comparing the larges't 

^5 of correlation results by magnitude witji a decision threshold. When Trexc^e'ds'thFpTesetT^^ 

Doppler frequency is taken as the estimate. Coherent aTidTlon-coherent processings are also incorporated in above- 
mentioned scheme to enhance the performance. The Doppler compensation prior to the correlation results In unde- 
graded correlation values when comparing with the techniques in which search for the Doppler frequency is performed 
after the correlation. The method according to the present invention is characterized in what will be presented in the 

50 characterizing part of the appended claim 1. The location system according to the present invention is characterized 
in what will be presented in the characterizing part of the appended claim 8. The receiver according to the present 
invention Is characterized in what will be presented In the characterizing part of the appended claim 15. The electronic 
device according to the present invention is characterized in what will be presented in the characterizing part of the 
appended claim 22. 

55 [0016] The present invention gives considerable advantages with respect to methods and receivers of prior art. The 
method of the invention can be used to improve and speed up the acquisition to the received signals also on weak 
signal conditions. By the method according to the invention, the acquisition calculation can be periormed by using less 
calculations than methods of prior art. The receiver according to the invention can be implemented with a relatively 
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20 


small number of components and the total energy consumption can be kept reasonable in comparison with a receiver 
of prior art to achieve the same resolution, wherein the invention is particularty well applbable in portable devices. 
Thus, the positioning receiver can also be implemented in connection with a wireless communication device. 
[00171 In the following, the present invention will be described in more detail with reference to the appended drawings, 
5 in which 

pjg la illustrates, in a reduced block chart, a signal monitoring block of a positioning receiver. In which 

* the method according to an advantageous embodiment of the invention can be applied, 

10 Fig. ib illustrates, in a reduced block chart, an acquisition block of a positioning receiver, according to 

an advantageous embodiment of the invention, 

Fjg 1 c illustrates, in a reduced block chart, the signal monitoring block of a positioning receiver, in which 

the method according to another advantageous embodiment of the invention can be applied, 

Pig Id illustrates, in a reduced block chart, an acquisition block of a positioning receiver, according to 

another advantageous embodiment of the invention, 

Figs. 2a and 2b illustrate the coarse and fine Doppler methods according to the invention. 

Fig. 3 illustrates the composition of received data with columns as repetitive epochs arranged in a row, 

Fig. 4 illustrates the correlation peak degradation as a function of the Doppler shift, 

Fig. 5 illustrates coherent addition in fine Doppler compensation phase according to an advantageous 

method of the invention with multiple assumptions on possible fine Doppler shift, 

Fig. 6 illustrates in-epoch compensation in fine Doppler compensation phase according to an advanta- 

geous method of the invention, 

Fig. 7 illustrates matched filtering in the fine Doppier compensation phase according to an advanta- 

geous method of the invention. 

Fig. 8 illustrates an advantageous embodiment of the non-coherent summing method which can be 

used in connection with the invention, 

Fig 9 illustrates a two-dimensional search matrix fomned by the non-coherent summing method which 

can be used in connection with the invention, 

Fig. 1 0 illustrates a two-dimensional search matrix formed by the shifting of the FFT transformed pseudo 

random sequence of the corresponding satellite according to another advantageous embodiment 
of the invention, and 

Fig. 1 1 illustrates an electronic device according to a preferred embodiment of the invention in a reduced 

block chart. 

[0018] The acquisition method according to an advantageous embodiment of the present invention will be described 
in more detail. The purpose of the acquisition is to adjust the phases of the incoming signal and the locally generated 
pseudo random sequence of the con^esponding satellite to a small timing offset and also find Doppler frequency mod- 
ulation. 

[0019] When the receiver is turned on and also in a situation in which the receiver has not been capable of receiving 
the signal of any satellite for a long time a two-dimensional search phase is performed in the receiver for each satellite 
whose signal is to be received. In this two-dimensional search phase the aim is to find out the shift in carrier frequency 
and the code phase of the satellites. The cannier frequency is affected by the Doppler shift and Inaccuracies of the local 
oscillator of the receiver. The uncertainty of the carrier frequency may be as large as ±6 kHz, wherein the receiver has 
to perform the search within 12 kHz frequency range around the transmission frequency (LI = 1575.42 MHz). Also, 
the receiver does not know the exact code phase wherein the receiver has to find out the con-ect code phase. The 
code period of a GPS signal consists of 1023 chips and when signal is sampled at a rate of k^ samples per chip then 
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there are ks'l 023 possible code phases for search. The number kg is not necessarily an integer, e.g.. 1 023 chips atter 
sampling could be represented by 1024 or 2048 samples. Then the two-dimensional search process is needed, in 
which one dimension Is the Doppler (frequency) shift within 12 kHz frequency range and the other dimension is the 
code phase from }^*^023 possible code phases. In the method according to the invention this two-dimensional search 
5 process is performed by using two phases: the coarse Doppler compensation phase and the fine Doppler compensation 
phase The coarse compensation phase utilises the whole frequency range, e.g. 12 kHz frequency range preferably 
in 1 kHz (or less) steps, and the fine Doppler compensation phase utilises the frequency range of each step of the 
coarse Doppler compensation, e.g. 1 kHz frequency range. It Is obvious that the values used in this disclosure are 
presented as non-limiting examples only. The present invention can also be applied with other than GPS systems. 
10 wherein the mentioned frequency values, code phases and number of codes may vary. 

[00201 In the following an advantageous method of the invention will be described m more detail tn connection with 
the receiver of Fig. 1a and the acquisition block of Fig. lb. To perfomn the coarse Doppler compensation phase along 
with the fine Doppler compensation phase the received signal is first sampled to convert the signal into discrete, co- 
herent values, i.e. complex numbers comprising digitized values of the Q-component and the l^omponent of the re- 
15 ceived signal. As was mentioned above the coarse Doppler search is performed in broader steps, e.g. irrl kHz steps 
and the fine Doppler search is performed in nan-ower steps within each coarse Doppler search, eio^ in 0.2 kHz steps. 
[00211 Fig. la presents schematically the acquisition process. The processing in acquisition b!o>>k.6^couIdJ^^ 
fonned at intermediate frequency with appropriate implementation of numerically controlled oscillator, but m the fol- 
lowing description it is assumed that the processing is perfomned at baseband. 
2o [00221 The frequency of the numerically controlled oscillator 4 is set such that the receiver receives the demanded 
frequencies of the frequency range to be examined. In this example the examined bandwidth is 12 kHz. Then the 
received signal is transfomied in a RF-block 2 into the intermediate frequency range or into baseband and sampled in 
an A/Dconverter block 3. The sampling rate of the A/D-converter block 3 is such that at least one coherent sample of 
every chip is taken. GPS C/A signal is periodic with the period of 1ms. In GPS system the chip rate is 1.023 chips (1 
25 epoch) in 1 ms wherein the sampling rate in acquisition stage is preferably 1024 samples in 1 ms, or an integer multiple 
of it Then a total of 1024 samples (or integer multiple of 1024 samples) are taken and stored into memory 5. For 
presentation purposes we will next assume that the number of samples is 1 024 .These samples form a sample vector 
X The 1024 samples of the sample vector represent a signal having ca. 1 ms length. The sampling of the signal 
continues wherein the next 1 024 samples are stored as the next sample vector x„. These sample vectors are continuous 
30 in time such that the following sample vector continues after the previous sample vector, i.e. the time difference between 
the last sample of the previous sample vector and the first sample of the next sample vector is substantially the same 
as the time difference between samples in the sample vectors. 

[0023] A number of sample vectors are stored for fine Doppler search. Depending on the granularity of the fine 
Doppler search and number of coherent additions the number of sample vectors may vary. For example, five consec- 

35 utive sample vectors are used wrth the fine Doppler search, which is illustrated as a matrix X in Fig. 3. The number of 
columns, i.e. the number of sample vectors x„ in the matrix X is labelled asUj and the number of rows, i.e. the number 
of samples x„ in each vector is labelled as N2 wherein the total number of elements in the matrix X is N-N^Ng. Although 
there are 1 024 samples indicated in Fig. 3 it is obvious that the number of samples can also be other than that. In one 
alternative sampling method 1023 samples (or multiple of 1023) are taken and then one (or multiple) extra values. 

40 preferably zeros, are added to the sample vectors to achieve a vector length which is a power of two. This makes the 
FFT algorithms easier compared with the situation where the length of the sample vectors is other than a power of two. 
[0024J A first time-to-frequency transfonn is performed after the sampling phase in the method according to the 
invention The first time-to-frequency transform is performed to the samples of the sample vectors such that in each 
transform one sample of each sample vector is used. For example, the first samples of each sample vector are used 

45 in one transfomi, the second samples of each sample vector are used in another transform, etc. Therefore 1 024 time- 
to-frequency transforms are performed for each fine Doppler search in this advantageous embodiment of the invention . 
The results of these time-to-frequency transfonns are saved as a first transfonn matrix Y. This time-to-frequency trans- 
form is preferably a fast Fourier transform (FFT) and it is illustrated as a first FFT block 12 in Fig. lb. 
[0025] Hereinbelow, the Fourier transfomi will be primarily used as an example of a time-to-frequency transform and 

50 an inverse Fourier transform as an example of an inverse transform, Le. a frequency-to-time transfomn; however, it is 
obvious that the present invention is not limited solely to these examples. 

[0026] An advantageous embodiment of the invention comprises also an in-epoch compensation phase 1 3 which Is 
performed after the first time-to-frequency transfomn phase to eliminate possible degradation of correlation peak values 
due to the Doppler shift. This degradation is illustrated in Fig. 4. A cun/e A illustrates the non-compensated correlation 
55 result in frequency domain. A cun/e B illustrates the compensated correlation result and a curve C illustrates the nor- 
malized compensated con-elation result when a windowing is used with FFT, respectively. From Fig. 4 it can cleariy be 
seen how the in-epoch compensation increases the accuracy of the correlation phase. This is achieved by multiplying 
14 the values of the matrix Y with respective values of a compensation matrix C (Fig. 6). The result of the potnt-by- 
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poinl muttiplication is illustrated as a compensated matrix CY in Fig. 6. 

[00271 result of the previous stage Is the fine Doppler compensated and coherently enhanced signal, meaning 

that each column of the matrix Y or CY is the Doppler compensated and coherently combined epochs of original signal 
assuming a particular .fine Doppler frequency. 

5 [0028] l4exX, the correlation phase is to be perfomied for each column of matrix Y or CY in which the correlation 
between the received signal of a satellite and the locally generated pseudo random code r of the satellite is calculated. 
The conelation phase is illustrated as a simplified diagram in Fig. 7a. The correlation result gives infomiation of the 
code phase of the received signal. The'correlation could be performed using methods of prior art but in the preferred 

I embodiment of the invention it is performed in frequency domain. 

10 [0029) The correlation in lime domain corresponds to a multiplication in frequency domain. Therefore, the matrix Y, 
^ or the compensated matrix CY. is transfomned 16 to frequency domain. The transfonn is performed in orthogonal 

direction with respect to the first time-to-frequency transform which in this advantageous embodiment of the invention 

means that the transform is performed by columns (epoch-by-epoch). 

[0030] The pseudo random codes r of each satellite of the. positioning system are preferably stored in the receiver, 
15 or the codes are generated in the receiver when needed. In the correlation phase the receiver selects, i.e. retrieves 
from the memory or generates the pseudo random code r which corresponds to the modulating code of the signal of 
the satellite into which the receiver is trying to perform the acquisition. This pseudo random code is called a replica 
code later in this description as a distinction to the actual pseudo random code used in the transmitting satellite. Then 
a time reversal is performed to the period of the replica code r. In this description the time reversal means such a 
20 modification that the first chip of the replica code becomes the last chip, the second chip becomes the second last 
chip, etc. The time reversed replica code is then transformed preferably to frequency domain by using e.g. fast Fourier 
transform. It is obvious that in some applications the time reversed version of the replica code can be stored into 
memory instead of the actual replica code r, or the receiver can directly generate the time reversed version of the 
replica code, wherein the time reversal of the replica code is not necessary during the operation of the receiver. It is 
25 also possible that if the time domain version of the reversed replica code is not needed, the transform domain version 
R of the reversed replica code can be stored in the receiver instead of the reversed replica code and the time-to- 
frequency transform of the reversed replica code is not needed during the acquisition phase. 

[0031] The matrix Y, or the compensated matrix CY, is directed as a first input into a correlator block 35 and the 
frequency domain version R of the reversed replica code is directed as a second input into the multiplier block 17. 

30 Then these two frequency domain inputs are multiplied and the result of the multiplication is transformed by columns 
into time domain in the frequency-to-time transform block 21 . It should be noted here that the matrices Y, CY, Zj are 
shown transposed In Fig. 7a, 8 for visualization purposes. The time domain output of the correlator block is then saved 
as a first correlation matrix (Fig. 8) for further inspection. Each column of the first correlation matrix Z^ comprises 
the crosscorrelation of the replica code r with all possible circular shifts of the same row of matrix Y, or CY. When the 

35 next correlation matrix Zj is formed in the similar way, it is added non-coherentty to the first one and the process 
continues for several consecutive matrices. 

[0032] Figure 8 illustrates the non-coherent addition of the correlation matrices Z^ — to form a search matrix S. It 
is evident, that the addition could be performed on-line without accumulating all the matrices. The non-coherent addition 
can be performed e.g. such that the absolute values of the elements of each correlation matrix Z are calculated and 
added up. The code phase and Doppler shift can now be searched by analyzing the search matrix S. In the analysis 
the values of the elements of the search matrix S are evaluated and advantageously a maximum value is searched. 
Then (he maximum value is found and if it exceeds the decision threshold the matrix indexes {row, column) of the 
maximum value define the correct code phase and Doppler shift which then can be used during the tracking and position 
calculation phases known as such. This non-coherent processing is illustrated as a block 22 in Fig. 1b. If the maximum 
^5 value does not exceed the decision threshold then it is decided that a signal is not present at this set of frequencies 
and the acquisition process continues with next coarse search stage. It is obvious that the above described method 
tor detecting a signal of a satellite is used as an example here but other types of detectors could be used for sign-ai 
detection as v/ell. 

[0033] The above described two-phase analyzing procedure is also illustrated in Fig. 2 where the lines on the left 
50 illustrate the frequency range division for the coarse Doppler search and the lines on the right illustrate the frequency 
range division for the fine Doppler search. 

[0034J The mathematk^l details of the invention are considered in the following. The sequence of received GPS 
data is denoted as S„ in a situation that the Doppler shift is not affecting to the signal. The number of samples in epoch 
is denoted as N2 and the number of periods during coherent processing is denoted as N^. Also some data addressing 
55 variables are defined: n^ , which indicates the epoch and n2, which indicate the position within the epoch. These variables 
are also illustrated in Fig. 3. The Doppler shift introduced by the satellite motion and local clock inaccuracies Is defined 
as k/(NiN2), i.e. taking the discrete grid of frequencies. Due to this Doppler shift the input signal is effected by a mod- 
ulation and an additive noise x^,-*^ in a form 
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[00351 Next, the multiplier ^ is neglected as it is a common scaling factor which it is not relevant for the further 
analysts and it cancels out during the non-coherent processing. In general. GPS receiver multiplies the incoming signal 
with a local replica code and integrates over several periods of the GPS signal with overall number of samples involved 

[00361 The coherent addition of epochs as the method to increase the signal-to-noise ratio will cause correlation 
peaks to be degraded if not compensated appropriately before the correlation as was mentioned earlier in the descrip- 
tion. The compensation can be performed by multiplying the incoming sequence by 

[00371 if the Doppler shift is already known. However, in practical situations the Doppler shift is normally not known 
a priori. One advantageous solution is to perform the compensation with different values of k wherein the correct 
frequency value will exhibit itself as an undegraded peak after the correlation. This is illustrated in Fig.2b. in which the 
method is illustrated as a simplified diagram. The input signal is examined in coarse frequency steps. The blocks D1, 
02, .... DM illustrate those coarse frequency steps. The coarse compensation is performed for each selected frequency 
step by mixing 201 the input signal with the signal of the numerically controlled oscillator 4. The frequency of the 
numericalfy controlled oscillator 4 is tuned 202 e.g. in 1 kHz steps for each examining phase. Then, the compensation 
for fine frequencies is performed 203 prior the correlation with replica code for different code phases 204. Then the 
correlation results are combined coherently 205 and non-coherently 206 to produce a search matrbc S. The search 
matrix X is examined 207 e.g. by comparing it with a preset threshold value. If the decision signal has at least one 
value which is greater than the threshold value, it can be assumed that there exists a signal of a satellite at the examined 
frequency range. Also other prior art detection methods can be used to decide if there exists a signal of a satellite at 
the examined frequency range. 

[0038] The variables k and n can be divided into two components such that k=k-,+k2Ni. and n=n^N2+n2, where k^, 
= 0,1,....A/^ -1, and Ac2,n2 = 0,r...,yV2''' ■ ^^"^^ ^ variable z is defined such that the components of replica vector 
r are involved: 


V -I -1 ' j\ 

Z= Z S ^n,V,+n,^ y ^ ^ J (3) 


[0039] If Doppler frequency quess is correct then 


55 


n,=0 * ' ^ n,=0 * n,=0 ' * ' n,=0 * rt,=0 
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where the periodicity of the signal x^Aj^+ng = is taken into account. If the replica code and incoming signal are 
aligned. I.e., = x^/ A = ±1, then 


/»js=0 n,-0 


[0040] It can be seen from the equation above that system gain is increased because of performing the addition of 
10 ' signals coherently, while the noise bandwidth is decreasing by the factor A/^Ng. The signal-to-noise ratio is also in- 
creasing by the same factor. 

[0041] An efficient method for calculating the sum 


/ij =0 n, «0 

20 in the equation (3) is described in the following. First, the inner sum is to be considered. From the splitting of k it can 
be seen that k, can be interpreted as the fine Doppler frequency shift while the k2 is the coarse Doppler frequency 
component. 

[0042] For each value of n2 and the frequency parameter k^ the coherent addition along will be 


25 


Af I .2/r, . 27r . / \ . 2/r . ■ 2;r, 

V - V r e ' ' e =e ^'^^ ■ v5 e (4) 

n.-O n,=0 


30 

where 


[0043] In the above equation for different values of k^ could be computed first by performing a time-to-frequency 
transfonn, for example discrete Fourier transform or fast Fourier transform, on rows of the matrix = f^^/MLg+n^ 
^0 then weighting the resulting matrix in a proper manner. Thus 


Y-=X (5) 

45 

where F is the time-to-frequency matrix and 

[0044] As can be seen from the equation (6) above, after the time-to-frequency transform of the rows, the elements 
of the matrix Y should be point-by-point multiplied with an appropriate complex sinusoid value. 
55 [0045] The fomnula 
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means thai the incoming signal x„ would be multiplied with a sinusoid defined by coarse frequency variable /Cj prior to 
the lime-to-frequency transfonn. This could be done before fine Doppler compensation as it is shown in Fig. la, Fig. 
lb and Fig. 2b but il could also be performed in a way presented in Figs. 1c and 1d without using of numerically 
controlled oscillator 4 for acquisition process. 
10 [00461 For the later case it is recalled that the meaning of the cross-correlation is to find the match in code phase 
between the received data and the replica code. Instead of compensating the coarse Doppler shift on the signal the 
replica code can be modulated by the same Doppler frequency: 


From the properties of the discrete Fourier transform the given modulation in time domain exhibits itself as a circular 
shift in frequency domain. Thus the frequency domain representation of the replica code can be shifted in a circular 
way without actual transformation for each coarse Doppler phase. The replica shifting approach is considered in detail 
later in this description. 

[0047) Replacing n with n^Nj+nj, x.„ can be represented as 


30 [0048] The equation (7) shows that the coarse acquisition modulation is independent on both the fine frequency 

and particular epoch n^. This means that the same modulation applies to all the epochs and it can be taken out of the 
time-to-frequency transfonn of the rows of the matrix X. In the following matrixes are used to further explain the cross- 
correlation. 

[0049J Next, the outer sum in the equation (3) is to be considered. In fact, this means actually that the replica code 
35 is correlated with the coherently enhanced incoming signal. This sum should be performed for all possible circular 
replica code shifts because the code phase of the incoming signal is unknown. Although other methods of prior art 
could be used the preferred embodiment Is an efficient approach to perform the cross-correlation for all possible shifts 
of replica code based on discrete Fourier transform {OFT) and its fast implementation (FFT). 

[0050] The cross -correlation method using DFT assumes that DFT of received signal fragment is point-by-point 
40 multiplied by the DFT of the time-reversed replica code and the inverse DFT shows the result of the cross-correlation. 
The following matrix formulation of this process can be used with a situation in which no Doppler compensation is 
performed: 

^5 Z^r' R FX (8) 

where Z is the cross-correlation matrix, in which each column I is the result of cross-correlation of the l-th column of 
X with the replica code, square matrix R is diagonal and diag{Fi) = F-'r^^ where r^^is the replica code vector. To include 
also the Doppler compensation the above equation (8) can be modified by point-by-point multiplication of two matrices 
C andx • F^. This point-by-point multiplication of two matrices is denoted as (.*) in this description. Then the cross- 
correlation can be written as 

55 Z= F'^ F-C(.*)(X F^) (9) 
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where matrix C is defined as 


30 


n2=0, 1 N2-I - The specified range of is due to the fact that the range of frequencies less than 1 kHz are to be 
considered at this stage. 

10051] When the coarse Ooppler frequency component described by the variable ka is taken into account, the equa- 
te tion (9) can be modified into a fonrt 

Z=F-^RF- C(.*)(m, ■XF^) = F-^RFMy- C{.*)(x - F^) (10) 

15 

where is a diagonal matrix and 

20 , , .X -;Tr~^2"2 


'N2 


So each column Is modulated arid then time-to-frequency transfonn such as FFT is perfonned. Taking into account 
25 that N2 is the epoch length and the properties of the DFT, modulation on the input by 

e 

circularly shifts the output of the time-to-frequency transfonn by ka. /.e- by an integer shift. If also the transform domain 
representation of the time-reversed replica code is circularly shifted, the result is a shifted vector prior to frequency- 
to-time transform, such as inverse FFT. Shifting in the transform domain can be considered as modulation in time 
domain. Therefore the coarse modulation on the input of the correlator exhibits itself as a modulation on its output if 
the transform domain representation of the time-reversed replica code is circularly shifted by kg. More fomnally, the 
cross-con-elation can be written as 


Z=M,-F-'-R,^,fl-F-C{.*)[x-F^) (11) 

Where R^hm 's the diagonal matrix with diagonal elements representing transform domain replica code shifted circularly 
by kg positions. The complex exponentials in are constant for each row and normally they are not signiHcant for 
non-coherent processing and detection. If necessary, they can be compensated by appropriate multiplication. As a 

45 conclusion from the mathematca! inspection above is that the modulation in the input (right of F in equation 9) is moved 
to the output (left of F in the equation 1 1 ). Therefore the first two time-to-frequency transforms in equation (1 1 ) are not 
necessary for each coarse frequency search. Fig.ld shows the structure of the acquisition block with the method of 
shifting the replica code. The numerically controlled oscillator 4 is not required for acquisition in this advantageous 
embodiment of the invention but it is used for tracking purposes. Another remarkable advantage of using the method 

50 of shifting the replica code is the possibility of significantly reducing the computation complexity at the expense of 
memory. In the method according to an advantageous embodiment of the invenUon the input signal is time-to-frequency 
transformed twice, in-epoch compensated, and stored into memory. This part is perfonned by the section 11 in Fig.ld. 
Then the result can be used for all coarse frequency stages and with all the replica codes. This will be performed by 
block 16 in Fig.ld. The complexity then is reduced almost two times and is mainly due to inverse DFT (FFT). 

55 [0052] Figure 10 illustrates an example of the two-dimensional search matrix S fomied by the shifting of the FFT 
transformed pseudo random repiica sequence of the corresponding satellite according to an advantageous embodi- 
ment of the invention. 

[0053] The operation of the acquisition block of the positioning receiver has only been presented above, but practical 
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receivers also comprise e,g. positioning functions which can comply with prior art as such. 

[0054] To calculate the position, the rece'rver performs the reception of a signal preferably on the basis of the signal 
received from at least four satellites. Thus, the above presented steps are repeated for the signal of each satellite, if 
necessary. 

[0055] A majority of the blocks required for applying the method can be implemented e.g. in a digital signal processor 
(not shown). To perform the FFT transforms, it is possible to use either hardware-based solutions or software applica- 
tions in the digital signal processor. Furthermore, to control the operation of the receiver it is possible to use a control 
means, preferably a microprocessor or the like, which is, however, not presented in the appended drawings, and which 
is prior art known to anyone skilled in the art. 

[0056] At least some of the blocks required for applying the method can also be implemented In a communication 
network (not shown), wherefore such blocks are not necessarily needed in the receiver In this kind of arrangement 
the receiver and the network are arranged to communicate with each other so that the steps according to the invention 
can be performed. 

[0057] The circular shifting of the replica code is also applicable in other FFT based correlation systems, where a 
phase shift of a received, code-phase modulated signal is to be detemnined, e.g. CDMA systems (code division-multiple 
access). 

[0058] Further, the appended Fig. 11 shows an electronk; device 23 complying with an advantageous embodiment 
of the invention, comprising functions of a wireless communication device and a positioning receiver. A first antenna 
24 is used to receive a signal transmitted from positioning satellites. The received signal is transferred to a first radio 
part 25, in which the signal is converted to an intemnediate frequency and digitized. The first radio part comprises e. 
g. a downconverter block 2, and an A/D-converter block 3 and a multiplier block 36 in the receiver 1 of Fig. la or Fig. 
1c. The digitized signal, which at this stage preferably comprises I and Q components, is transferred to a digital signal 
processing unit 26, in which e.g. the above-presented correlations to the different phases of the replica code signal r 
are perfomned. The electronic device 23 also comprises first memory means 27 which comprise e,g, a random access 
memor/ for storing data required during the operation of the receiver as well as preferably also a read only memory 
and/or an non-volatile random access memory for storing the program code of the digital signal processing unit 26. In 
this embodiment, for example part of the functions of the control block 7, such as the FFV transform blocks 12, 15. the 
multipliers 14, 17, 20 and the inverse FFT block 21 , are implemented in the signal processing unit 26. The digital signal 
processing unit 26 comprises preferably also means for generating and shifting the frequency domain replica code. It 
is obvious that at least some of said functional blocks can also be implemented in another way known per se. For 
example, the FFT transform blocks 12, 15 and the inverse FFT transform block 21 can also be implemented with 
integrated circuits, or the like, designed for these purposes. 

[0059] The digital signal processing unit 25 transfers the values for determining the code phase and the frequency 
deviation to be used in the determining block 22 to a processor block 28 which comprises e.g. a microprocessor and 
I/O logics. Some of the functions of the acquisition block 6 are preferably at least partly implemented as program 
commands in the processor block. As the data memory and the program memory for the processor block 28, second 
memory means 29 arc used. It is obvious that the first memory means 27 and the second memory means 29 can also 
comprise a common memory. Positioning information can be displayed to the user on a display 30. 
[0060] Also functions of the wireless communication device are implemented in the application software of the proc- 
essor block 28. Thus, the display 30 can also be used, in a way known perse, to display e,g, call information. Using 
a keypad 31 , the user can control the positioning receiver as welt as the wireless communication device. A codec 32 
is used to encode and decode audio signals. Further, Fig. 11 shows a radio part 33 and a second antenna 34 for the 
wireless communication device. 

[0061] The present invention is not limited solely to the above presented embodiments, but it can be modified within 
the scope of the appended claims. 

Claims 

1- A method for acquiring a receiver (1) into a code modulated spread spectrum signal received by the receiver (1), 
in which method at least one replica code (r) is used which corresponds to a code used with the modulation having 
a pre-determined number of chips, and an examination phase is performed. In which a frequency shift of the signal 
wilhin a selected frequency area is examined, and a code phase of the code used with the modulation is examined, 
characterized in that the frequency shift examination is divided into a first estimation phase and a second esti- 
mation phase, wherein in the first estimation phase the selected frequency area is divided into a first set of fre- 
quencies, and in the second estimation phase a second set of frequencies is examined nearby each frequency of 
the first set of frequencies, that a comparison for frequencies of the second set of frequencies is performed using 
the received signal and the replica code (r), and that the results of the comparison are used to estimate the correct 
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frequency shift. 

2. The method according to the claim 1 , in which a reference oscillator signal is formed, characterized In that in 
each of the first estimation phases the received signal (xj is mixed with said reference oscillator signal, that the 

5 frequency of the reference oscillator is set into a different frequency for different first estimation phases, and that 

the mixed signal is used in the second estimation phase. 

3. The method according to the claim 1 , characterized In that a time-to-frequency transformation of a reversal of 
the replica code (r) is produced, that in each of the first estimation phase the transfomned, reversed replica code 

10 (r) is shifted such that in different first estimation phases a different phase shift of the transfomned, reversed replica 

code is used- 


4. The method according to any one of the claims 1 , 2 or 3, characterized In that the received signal is sampled for 
producing a set of samples (x^), a matrix (X) is formed from the samples, the matrix (X) having a first dimension 
and a second dimension, which first dimension preferably equals the number of samples of the code period, and 
that the second estimation phase having the steps of performing a first time-to-frequency transform' on the matrix 
(X) in said second direction, and performing a second time-to-frequency transform on the time-to-frequency trans- 
formed matrix (X) in said first direction. 


20 5. The method according to any one of the claims 1 , 2 or 3, characterized in that the received signal is sampled for 
producing a set of samples (xj, a matrix (X) is formed from the samples, the matrix (X) having a first dimension 
and a second dimension, which first dimension preferably equals the number of chips of the code, a compensation 
matrix (C) is fomied, and that the second estimation phase having the steps of performing a first time-to-frequency 
transform on the matrix (X) in said second direction, multiplying the time-to-frequency transformed matrix (X) with 

25 the compensation matrix (C) to fomri a compensated matrix (CX), and performing a second time-to-frequency 

transform on the compensated matrix (CX) in said first direction. 

6. The method according to any one of the claims 4 or 5, characterized in that a time-to-frequency transform is 
perfomned on the reversed replica code (r), the time-to-frequency transfonmed replica code (R) is multiplied with 

30 the resulting matrix of the second time-to-frequency transformation, a frequency-to-time transfomi is performed 

on the resulting matrix of the multiplication. 

7. The method according to the claim 6, characterized in that a non-coherent processing is performed on at least 
one frequency-to-time transformed matrix, in which non-coherent processing a maximum value is searched for 

35 finding a correct frequency shift and code phase. 

8. A location system comprising at least 

a receiver (1 ) having means (24, 25) for receiving code modulated spread spectrum signal, 

- means (6) for acquiring the receiver (1 ) into the received signal, 

means (1 8) for using at least one replica code (r) which corresponds to a code used with the modulation, which 
code having a pre-determined number of chips, and 

examination means (4. 6, 36) for examining a frequency shift of the signal within a selected frequency area, 
and a code phase of the code used with the modulation, 

45 

Characterized in that the examination of the frequency shift is divided into a first estimation phase and a second 
estimation phase, wherein the location system further comprises 

- means (4, 18, 36) for dividing the selected frequency area into a first set of frequencies for the first estimation 
phase, 

means (1 1 ) for examining a second set of frequencies nearby each frequency of the first set of frequencies in 
the second estimation phase, and 

means (16) for performing a comparison for frequencies of the second set of frequencies by using the received 
signal and the replica code (r), 

55 

and that the examination means comprise means (21 , 22) for evaluating the correct frequency shift by using the 
results of the comparison. 
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9. The location system according to the claim 8 further comprising a reference oscillator (4) for producing a reference 
oscillator signal, characterized In that the location system further comprises 

means (36) for mixing the received signal (x^) with said reference oscillator signal in the first estimation phases, 
5 and 

means (7. 8, SW2) for adjusting the frequency of the reference oscillator into a different frequency for different 
first estimation phases, 

and that the mixed signal is arranged to be used in the second estimation phase. 

10 

10. The location system according to the claim 8, characterized In that it comprises means (1 8) tor producing a time- 
to-frequency transformation of a reversal of the replica code (r), and for shifting the transformed, reversed replica 
code (r) in each of the first estimation phase such that in different first estimation phases a different phase shift of 
the transformed, reversed replica code is arranged to be used. 

15 

11. The location system according to any one of the claims 8, 9 or 10, characterized in that it comprises 

means (5) for sampling the received signal for producing a set of samples (x^), 

means (6) for forming a matrix (X) from the samples, the matrix (X) having a first dimension and a second 
?o dimension, which first dimension preferably equals the number of samples of the code period, 

means (12) for performing a first time-to-frequency transform on the matrix (X) in said second direction, and 
means (1 5) for performing a second ttme-to-frequency transform on the time-to-frequency transformed matrix 
(X) in said first direction. 


25 12- The location system according to any one of the claims 8. 9 or 10. characterized in that it comprises 
means (5) for sampling the received signal for producing a set of samples {x„), 

means (6) for forming a matrix (X) from the samples, the matrix (X) having a first dimension and a second 
dimension, which first dimension preferably equals the number of samples of the code period, 
^0 - means (13) for forming a compensation matrix (C), 

means (14) for multiplying the time-to-frequency transformed matrix (X) with the compensation matrix (C) to 
form a compensated matrix (CX), and 

means (1 5) for perfonning a second time-to-frequency transfonn on the compensated matrix (CX) in said first 
direction. 

35 

13. The location system according to any one of the claims 11 or 12, characterized in that it comprises 

means (18) for forming a time-lo-frequency transformed reversed replica code (R). 

means (14) for multiplying the time-to-frequency transformed reversed replica code (R) with the resulting matrix 
of the second time-to-frequency transformation, and 

means (15) for performing a frequency-to-time transform is perfonned on the resulting matrix of the multipli- 
cation. 

14. The location system according to the claim 13, characterized In that it comprises means (18) for performing a 
^5 non-coherent processing on the frequency-to-time transformed matrix, in which non-coherent processing a max- 
imum value is arranged to be searched for finding a correct frequency shift and code phase. 

15. A receiver (1) comprising at least 


means (24. 25) for receiving code modulated spread spectrum signal, 
means (6) for acquiring the receiver (1) into the received signal, 

means (1 8) for using at least one replica code (r) which corresponds to a code used with the modulation, which 
code having a pre-determined number of chips, and 

examination means (4, 6. 36) for examining a frequency shift of the signal within a selected frequency area, 
and a code phase of the code used with the modulation, 

characterized In that the examination of the frequency shift is divided into a first estimation phase and a second 
estimation phase, wherein the receiver (1) further comprises 
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- means (4.18, 36) for dividing the selected frequency area into a first set of frequencies for the first estimation 
phase, 

- means (11) for examining a second set of frequencies nearby each frequency of the first set of frequencies in 
the second estimation phase, and 

5 - means ( 1 6) for performing a comparison for frequencies of the second set of frequencies by using the received 

srgnal and the repiica code (r), and that the examination means comprise means (21, 22) for evaluating the 
correct frequency shift by usmg the results of the comparison. 

16. The receiver (1) according to the claim 15 further comprising a reference oscillator (4) for producing a reference 
10 oscillator signal, characterized In that the receiver (1) further comprises 

- means (36) for mixing the received signal (x„) with said reference oscillator signal in the first estimation phases 
and 

- means (7. 8, SW2) for adjusting the frequency of the reference oscillator into a different frequency for different 
first estimation phases, 

and that the mixed signal is arranged to be used in the second estimation phase. 

17. The receiver (1) according to the claim 15. characterized in that it comprises means (18) for producing a time- 
to-frequency transformation of a reversal of the replica code (r), and for shifting the transformed, reversed replica 
code (r) in each of the first estimation phase such that in different first estimation phases a different phase shift of 
the transformed, reversed replica code is arranged to be used. 

18. The receiver (1) according to any one of the claims 15. 16 or 1 7, characterized in that it comprises 

- means (5) for sampling the received signal for producing a set of samples (xj. 

- means (6) tor forming a matrix (X) from the samples, the matrix (X) having a first dimension and a second 
dimension, which first dimension preferably equals the number of samples of the code period 

- means (12) for perfonming a first time-to-frequency transform on the matrix (X) in said second direction and 

- nrieans (15) for perlomiing a second time-to-frequency transform on the time-to-frequency transformed matrix 
(X) in said first direction. 

19. The receiver (1) according to any one of the claims 15. 16 or 17. characterized in that it comprises 

- means (5) for sampling the received signal for producing a set of samples (x ) 

' means (6) for forming a matrix (X) from the samples, the matrix (X) having a'first dimension and a second 
dimension, which first dimension preferably equals the number of samples of the code period 
means (13) for forming a compensation matrix (C), 

- means (14) for multiplying the time-to-frequency transformed matrix (X) with the compensation matrix (C) to 
" 'ofrn a compensated matrix (CX), and 

" direction ^""^ Pertomr^ing a second time-to-frequency transform on the compensated matrix (CX) in said Hrst 

20. The receiver (1 ) according to any one of the claims 1 8 or 1 9. characterized fn that it comprises 

- means (18) for forming a time-to-frequency transformed reversed replica code (f^) 

- nieans (14) for multiplying the time-to-frequency transformed reversed replica code (R) with the resulting matrix 
01 the second ttme-to-frequency transformation, and 

^0 ' caiZ^ performing a frequency-to-tlme transform is performed on the resulting matrix of the multipli- 

21. The receiver (1) according to the claim 20. characterized In that it comprises means (18) for performing a non- 
coherent processing on the frequency-to-time transformed matrix, in which non-coherent processingTmaximum 
value IS arranged to be searched for finding a correct frequency shift and code phase. maximum 

22. An electronic device (23) comprising at least 

- a receiver (1 ) having means (24. 25) for receiving code modulated spread spectrum signal. 
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. means (6) for acquiring ihe receiver (1 ) into the received signal, 

- means (18) for using at least one replica code (r) which corresponds to a code used with the modulation which 
code having a pre-cJeiermined number of chips, and 

- examination means (4. 6. 36) for examining a frequency shift of the signal within a selected frequency area 
5 and a code phase of the code used with the modulation. 

characterized In that the examination of the frequency shift is divided into a first estimation phase and a second 
estimation phase, wherein the electronic device (23) further comprises 

' phase ' ^' '''' ""^''''''"^ frequency area into a first set of frequencies for the first estimation 

- means (11 ) for examining a second set of frequencies nearby each frequency of the first set of frequencies in 
the second estimation phase, and 

- means ( 1 6) for performing a comparison for freq uencies of the second set of frequencies by using the received 
signal and the replica code (r), 

and that the examination means comprise means (21. 22) for evaluating the correct frequency shift by using the 
results of the comparison. ' y 

20 23. The electronic device (23) according to the claim 22, characterized in that ft further comprises means (33 34) for 
communicating with a mobile communication network. \ . ; 
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